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Abstract 

We investigate the production of e + e~ pairs from p + Be and nucleus-nucleus 
collisions from 2 GeV/A to 200 GeV/A within a covariant transport approach and 
focus on the dilepton angular anisotropies as a function of the dilepton invariant 
mass. Whereas the low mass regime yields information about the Dalitz decays 
of the A, r\ and u, above M «0.5 GeV the angular anisotropies provide additional 
information about the it + it~ — > p° — > e + e~ channel in the medium. The anisotropy 
coefficient is found to change its sign for M > 0.5 GeV in case of nucleus- nucleus 
reactions when increasing the bombarding energy from 2 GeV/A to 200 GeV/A 
which sheds some light on the irir — * p dynamics. 
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1 Introduction 



Nowadays, dileptons are used as electromagnetic signals from the hot and dense nuclear 
phase in heavy-ion collisions at BEVALAC/SIS 0], @, | or SPS energies g § |. The 
information carried out by leptons may tell us not only about the interaction dynamics 
of colliding nuclei, but also on properties of hadrons in the nuclear environment or on a 
possible phase transition of hadrons into a quark-gluon plasma (cf. 0). However, there 
are a lot of hadronic sources for dileptons because the electromagnetic field couples to all 
charges and magnetic moments. In particular, in hadron-hadron collisions, the e + e~ pairs 
are created due to the electromagnetic decay of time-like virtual photons which can result 
from the bremsstrahlung process or from the decay of baryonic and mesonic resonances 
including the direct conversion of vector mesons into virtual photons in accordance with 
the vector dominance hypothesis. In the nuclear medium, the properties of these sources 
may be modified and it is thus very desirable to have experimental observables which 
allow to disentangle the various channels of dilepton production. 

Apart from the differential e + e~ spectra the investigation of lepton angular distribu- 
tions is promising, too, because the virtual photon created in hadronic interactions is 
polarized. The coupling of a virtual photon to hadrons induces a dynamical spin align- 
ment of both the resonances and the virtual photons due to the conservation laws and 
consequently the angular distribution of a lepton pair with respect to the polarized photon 
momentum is anisotropic. 

The angular distributions of leptons for large invariant masses (above the J/^> regime) 
were investigated at high energy, both experimentally || and theoretically (cf. the review 
of Strojnowski ||), more than a decade ago. Here, the dominant source at large invariant 
masses is the quark annihilation (Drell-Yan) process. The former measurements allowed 
to observe a small deviation of the polar angular distribution from the prediction of 
the naive quark-parton model, the deviations observed, however, are in line with more 
detailed QCD computations. Thus the sensitivity of lepton angular characteristics already 
has been successfully used to differentiate between models. 

More recently, it has been proposed |T0| , |TT| to use the lepton pair angular distributions 
also for a distinction between different sources in the "small" invariant mass region (M < 
1 GeV), where a lot of dilepton sources contribute. It has been shown, furthermore, that 
due to the spin alignment of the virtual photon and the spins of the colliding or decaying 
hadrons, this lepton decay anisotropy turns out to be sensitive to the specific hadronic 
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production channel. In our previous works [10[ [11], |12[ we have calculated the dilepton 
anisotropy for elementary nucleon-nucleon collisions and for p + d reactions at BEVAL AC 
energies. First results for proton-nucleus and nucleus- nucleus collisions at 1 - 2 GeV/A 
are reported in |H| , where the calculated inclusive dilepton spectra were compared to the 
DLS data 0, g |. 

In this work we present a more systematic study on the dilepton angular anisotropy 
for p + A and A + A reactions from 2 GeV/A to 200 GeV/A within the covariant transport 
approach HSD |L4|]. At SPS energies the inclusive dilepton spectra can be compared to the 
results of the CERES/NA45 and HELIOS-3 collaborations @, ||, a task that already 



has been addressed with some success in |L5|, [16[ . Here we proceed with predictions for the 
corresponding dilepton angular distributions and for the system Au+Au at 10 GeV/A and 
160 GeV/A. Our work is organized as follows: in Section 2 we describe the evaluation of 
the lepton pair angular distribution and briefly repeat the main ingredients of the covariant 
transport approach. Section 3 contains the calculated results on dilepton spectra for p+Be 
and various nucleus-nucleus collisions as well as the dilepton anisotropics as a function 
of the invariant dilepton mass. We summarize our study in Section 4 and present a brief 
discussion on the experimental requirements to measure the effects proposed. 



2 Theoretical ingredients 

2.1 Definition of the anisotropy coefficient 

The general form of the angular distribution for the decay of a virtual photon into a lepton 
pair may be written 



as 



3 r 

W(6, <p) = — Pn(l + cos 2 9) + (1 - 2p u ) sin 2 9 + pi_i sin 2 9 cos 2 2(p 



8tt 



+V / 2Repio sin 29 cos 



where the angles 9, if of the electron momentum are measured with respect to a fixed 
z-axis in the virtual photon rest frame. The density matrix elements pij depend on the 
choice of the reference frame as well as all the variables describing the virtual photon. An 
integration over the azimuthal angle gives 

W{9) ~ l + Bcos 2 #, (2) 
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where the coefficient B may vary between —1 and +1. Adopting the normalization con- 
dition poo + 2pn = 1, the coefficient B can be represented as (cf. [fill): 

B = (3) 

1 - Pll 

Adopting the general form Eq. ([!]) the angular distribution of dileptons created in 
channel i in a hadron-hadron (h + h) reaction then can be written as [fL0[ , |i~T|1 

da} h 



dMd cos 0, 



/i/i 



A^(l + B? h cos 2 9 hh ). (4) 



Here 0^ is defined by 9hh = (l-,v% h ) with the electron momentum l*_ measured in the 
dilepton center-of-mass system (g* = l*_ + 1+ — 0), while = ^ jq^ 1 is the velocity of 
the dilepton c.m.s. relative to the h + h cm. system. 

The coefficient B\ h in Eq. (|j) describes the anisotropy of the angular distribution 
(B^ = in case of isotropy), while A\ h determines the magnitude of the respective cross 
section. The absolute value of B\ now depends on the choice of the coordinate system. 
For example, in the Gottfried- Jackson frame the z-axis coincides with the direction of the 
incident beam in the dilepton rest frame; this definition was adopted in the experiments 
before ||. Following |T0| , |i~Tf we will use for our analysis the "helicity" system, where the 
z-axis is chosen along the direction of the virtual photon momentum q in the center-of- 
mass system (c.m.s.) of the colliding hadrons. 

The total differential cross section for h + h collisions now can be represented as a sum 
of the differential cross sections for all channels i, 

A Jf\ = E j J? n = A h \M){l + B h \M) C o^6 hh ), (5) 
dMd cos hh ■ u , dMd cos hh 

mi i=ctiannel 

which leads to the total anisotropy coefficient, 

daf B™ 1 



dM i + I rWi 

B h \M)= £ <B? h (M)>, <B{ h (M)>= - hh +3 / , (6) 



^channel 



y-, dn 



dM ill B hh 

x ~ 3 « 



where the special weighting factors originate from the necessary angle-integrations. Thus, 
the anisotropy coefficient B hh for h + h reactions is the sum of the "weighted" anisotropy 
coefficients (< B\ h >) for each channel i obtained by means of the convolution of B\ h 
with the corresponding invariant mass distribution (cf. |]12[). 
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The anisotropy coefficients for the bremsstrahlung, Dalitz decays of A-resonance and 
??-meson in pn and pp interactions have been calculated on the basis of a one-boson- 



exchange model fitted to elastic NN scattering in Ref. [11]. These anisotropy coefficients 



are a function of the invariant mass M, the masses m a , mi, and the initial invariant energy 



y/s of the hadrons a + b involved in the reaction. As shown in Ref. []T0 1 , the anisotropy 
coefficient for pion annihilation in the 7r + 7r~ c.m.s. is given by B^+ v - = —1. For this 
particular situation, where there are only two particles in the initial and final states, 8hh 
is the angle of the lepton momentum with respect to the pion momentum in the c.m.s. 
of the leptons (or pions), i.e. q* = p* a + p* h = l* + + l*_ = 0. 

2.2 Covariant transport approach for nucleus-nucleus collisions 

The dynamical evolution of proton-nucleus or nucleus-nucleus collisions is described by 
a coupled set of transport equations evolving phase-space distributions fh(x,p) for each 
hadron h [[14] 



i.e. 



{(n, - u u d^ h - M^c/f) as + fa&pK + M* h d*u£) %} f h (x, P ) 
= ]T |d2d3rf4...[G t G] 12 ^34...5r(n + n 2 -n 3 -n 4 ...) 

X {fh :i ( x ^P3)fh i (x,PA)f h (x,p)f h2 (x,p 2 ) 

- fh(x,p)fh 2 (x,p 2 )fh 3 (x,P3)f h4 (x,p 4 )} . . . . (7) 

In Eq. (|7[) Ufr(x,p) and U£(x,p) denote the real part of the scalar and vector hadron 
selfenergies, respectively, while [G + G]i2^34...5 4 (IT + n 2 — n 3 — n 4 . . .) is the 'transition 
rate' for the process 1 + 2— > 3 + 4 + . . . which is taken to be on-shell in the semiclassical 
limit adopted. The hadron quasi-particle properties in (]?]) are defined via the mass-shell 
constraint |T9|| , 

5{Yi^-Mf) , (8) 
with effective masses and momenta given by 

M*(x,p) = M h + U*(x,p) 

W(x,p) = j?-U£{x,p) , (9) 

while the phase-space factors 

f h (x,p) = l±f h (x,p) (10) 
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are responsible for fermion Pauli-blocking or Bose enhancement, respectively, depending 
on the type of hadron in the final/initial channel. The dots in Eq. (0) stand for fur- 
ther contributions to the collision term with more than two hadrons in the final/initial 
channels. The transport approach (0) is fully specified by U^(x,p) and U%(x,p) (// = 
0, 1,2,3), which determine the mean-field propagation of the hadrons, and by the tran- 
sition rates G^G 5 4 (. . .) in the collision term, that describe the scattering and hadron 
production/absorption rates. 

In the HSDQ approach the following baryon states are explicitly propagated: the nu- 
cleon, A, A/"* (1440), N* (1535), A, £ and 5 as well as their antiparticles with proper isospin. 
The meson sector includes all pseudo-scalar mesons as well as the 1~ vector- meson octett. 
For the baryon selfenergies we use the covariant parametrizations described in detail in 



14jl , whereas all mesons (m) in this work will be propagated without selfenergies, i.e. 



= = 0. The inelastic transition rates are as described in Ref. |14|], i.e. for invari- 
ant collision energies \/s < 2.6 GeV we employ explicite cross sections for baryon-baryon 
reactions as in the BUU approach of Ref. |20| whereas for y^s > 2.6 GeV we employ the 



LUND string formation and fragmentation model ]2TJ] to describe the hadron production 



cross sections. As in |T4| a string formation time r = 0.8 fm/c has been used to translate 
the elementary 'free' cross sections to hadron reaction rates. In case of meson - baryon 
reactions explicite cross sections - in line with the available experimental data - are taken 
into account for y/s < 1.8 GeV whereas the LUND model is used above 1.8 GeV invariant 
collision energy. Meson - meson reactions, furthermore, are treated as binary collisions 
including Breit-Wigner resonance cross sections with resonance properties taken from [^2j . 
The coupled set of transport equations (0) is solved within the testparticle method 



employing the parallel emsemble algorithm [33] , where a hadron is represented by a single 
testparticle in each individual parallel ensemble. In this way the 4-momenta of all hadrons 
are known throughout the collision together with their inelastic and elastic transitions. 
We note that the HSD approach quite well describes the reaction dynamics for p + A and 
A + A collisions from SIS to SPS energies as demonstrated in JTjj] . 

The production of dilepton pairs (e + e~) includes the Dalitz decays A — > Ne + e~,r] —>■ 
7e + e~, to — > 7r°e + e~, the bremsstrahlung for charged AW, irN and tcti collisions as well as 
the direct decays of the vector mesons p, uo and <fi. Furthermore, the secondary mesonic 
channels 7t + tt~ — > p° — > e + e~,KK — ► — > e + e~ and np — > — > e + e~ are taken into 



^^Hadron-String-Dynamics 
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account, too. The decay ratios are taken from the particle data table p2] whereas the 
formfactors for the Dalitz decays are adopted from Landsberg [24]. Further details of the 
computations are given in Refs. |2(J [15], [16] and don't have to be repeated here. 



2.3 Anisotropy coefficients for p + A and A + A reactions 

For heavy-ion reactions the situation becomes more complicated due to the nuclear dy- 
namics and the explicit time evolution of the interacting system. For the calculation 
of the nuclear anisotropy coefficient we start from the point that the form of the angu- 
lar distribution for all "elementary" hadron-hadron interactions a + b, that occur in the 
nucleus-nucleus reaction A + B, are known from the microscopic calculations |T0| , |i~T| , |i~2|l . 

In case of A + A reactions we have to take into account that the "elementary" a + b 
c.m.s. (in which the Bf b are computed) are moving relative to the A + B c.m.s. This 
implies that the direction of the virtual photon momentum q ab relative to the hadron- 
hadron c.m.s. a + b is different for each elementary interaction a + b and can't be restored 
from the experimental data. However, the direction of the virtual photon momentum q AB 
in the nucleus-nucleus c.m.s. A + B can be reconstructed for each lepton pair. Thus, it is 
necessary to perform an angular transformation from 6 a b to 8ab, where &ab is the angle 
between the lepton momentum I* and the velocity of the dilepton c.m.s. relative to the 
c.m.s. of the colliding nuclei A + B {v q = q AB /Qq B )- 

At high energies (~ 200 GeV/A) the main contributions to the dilepton spectra come 
from the Dalitz decays of rj, tu-mesons, 7r + 7r~-annihilation and direct decay of vector 



mesons p,u,(j). As shown in Ref. |L5| the contributions from the A-Dalitz decay and 
bremsstrahlung are negligible. Furthermore, the polarization of vector mesons created in 
direct nucleon-nucleon reactions is practically zero according to the measurements of Blo- 
bel et al. ||25|| . This is due to the fact that the exclusive cross section NN — * NNp°, u, <fi 
is very small compared to the inclusive cross section iViV — > NNp ,u,<j) + X containing 
several pions in the final channel, too. Consequently, we do not consider the contribu- 
tion to the anisotropy coefficient from the decays of vector mesons produced in primary 
iViV-collisions. Thus, there are only three dominant channels at intermediate and high 
energies, i.e. the Dalitz decays of rj - and to - mesons as well as 7r + 7r~ annihilation that 
have to be taken into account in our analysis. 

For the computation of the nuclear anisotropy coefficient we use another algorithm 



than in our previous paper [13|, i.e. we simulate explicitly lepton events for each channel 
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i with fixed angular distribution in the individual a + b system and then transform the 
lepton and photon 4-momenta to the A + B system. For example, for the Dalitz decay 
of the ?7-meson, the lepton angular distribution in the rest frame of the ?7-meson has a 
form @) with B ab = +1 [H|, where 9 a b is the angle between the lepton momentum I* 
and the direction of (f 7 , while <f is the virtual photon momentum in the rest frame of 
the rj. We also use that the distribution of the virtual photon momentum <f in the rest 
frame of the rj is isotropic. Thus, for each rj event with energy E v and momentum P v 
in the nucleus c.m.s. A + B we generate lepton events distributed with respect to <f 
according to Eq. (Q) with B^ b = +1 in the lepton c.m.s. The next step then is the Lorentz 
transformation from the r] rest frame to the c.m.s. of A + B: (Jab = L {^E~^j where 

1°v = "~2^r"~' \$n\ = \i ~2m' '°' > " Fi nan y we compute the angle Qab between I* and v AB 



as 



I* ■ v AB 

008043 = — q - , (11) 

\l*\.UAB\> 



where v AB = q AB /q$ B and |/*| = M/2 for fixed invariant mass M. Thus the angular 
distribution as a function of cos6Ub for all generated lepton events is recovered and B v 
can be computed via 

f ab _ W(9 AB = 0°) , , 

^ W(6 AB = 90°) ' 1 ; 

The weighted < B v > then follows from Eq. (f|). 

In a similar way we calculate the anisotropy coefficient for the pion annihilation channel 
by using that the anisotropy coefficient in the c.m.s. of the leptons (or pions) JK| B^ b +7T - = 
— 1. For each pion annihilation event we then restore the virtual photon energy and 
momentum in the c.m.s. of A + B: q AB = E a + Eb, q AB = p a +Pb, where E a , Eb, p a ,Pb are 
the energies and momenta of the pions in the A + B c.m.s. Then we perform a Lorentz 
transformation for the pion momentum from the A + B c.m.s. to the leptons c.m.s.: 
p* a = L {^XB^jPa- In the lepton c.m.s. we generate the lepton events with |/*| = M/2 
(M 2 = q As2 — q AB2 ) and angular distribution W(9 a b) = 1 — cos 2 6 a b, where Q a b is the angle 
between the I* and the pion momentum p* a in the leptons c.m.s.; cos 6 a b = J . . For 

r \'\Pa\ 

each selected lepton event the angle 9ab can be calculated according to Eq. flTT|). The 
weighted coefficient < B^^- > can be computed in the similar manner as in the rj case. 

The anisotropy coefficient for the Dalitz decay of the tu-meson, furthermore, is com- 
puted in analogy to the rj case, using the elementary coefficient B® b = +1. 
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We have tested our algorithm for BEVALAC energies and - within the statistical 
accuracy - reproduce the results from our previous paper where an approximation 
was made to calculate the anisotropy coefficients for the Dalitz decays of A and 77: in |L3 



we have used the results of our microscopic calculations for the anisotropy coefficient that 
were already averaged over q Clb whereas, in principle, we need this degree of freedom to 
transform correctly the differential cross section from the a + b to the A + B c.m.s. Our 
discrete algorithm now allows to avoid this approximation; it is close to an experimental 
event by event analysis and includes the full hadron dynamics. 



3 Numerical results for p + Be and A + A reactions 

In this section we present the results of our calculations for the dilepton spectra and the 
anisotropy coefficient for the systems p+Be, Ca+Ca, S+Au, Au+Au from BEVALAC to 
SPS energies without employing any medium effects for the mesons, i.e. = 0, = 
in the set of transport equations (0) for all mesons. 

In the Fig. [l] we show the computed weighted anisotropy coefficients < Bi(M) > for 
p+Be and Ca+Ca collisions at the bombarding energy of 2 GeV/A essentially confirming 
our previous results from |13| employing the averaged coefficients as described above. We 



note that the inclusive dilepton spectra for the reactions have been presented in Ref. JL3 
together with a "cocktail" decomposition. The main contributions arise from the r\ and 
A Dalitz decays due to their large 'elementary' anisotropy coefficients and cross sections, 
respectively. The contribution from pn bremsstrahlung is practically zero at all invariant 
masses due to a small 'elementary' anisotropy coefficient. The weighted coefficient from 
7r + 7r~ annihilation is rather small (~ 0.1) even for the Ca + Ca reaction and decreases 
for M > trip due to the threshold behaviour of the cross section. However, compared to 
< B n + n - (M) > for p + Be, where pion annihilation is very low, a clear (but moderate) 
enhancement can be extracted. 

Before going over to the presentation of the anisotropy coefficient at AGS energies 
(~ 10 GeV/A) we show the differential multiplicity dn e + e -/dM for p + Be and Au + Au 
in order to demonstrate the relative contribution of the different channels employing an 
'energy resolution' of 50 MeV (Fig. |^). In the 'cocktail' plot the u — > 7re + e~ Dalitz decay, 
rj Dalitz decay (V), the pion annihilation channel ('7r + 7r~') as well as the direct decay 
of the vector mesons ('p, u;, $') are displayed explicitly while the solid curves (denoted 
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by 'a//') represent the sum of all sources. The contributions from the A Dalitz decay 
and the bremsstrahlung channels are low by an order of magnitude compared to the 
?7-channel and not plotted explicitly. Similar to the systems at BEVALAC energies the 
7r + 7r - annihilation channel is of no importance for p+Be, but becomes the most significant 
sources for M > 0.4 GeV in case of Au + Au collisions. Furthermore, when scaling the 
spectra for Au + Au by the number of projectile nucleons, we observe an additional 
enhancement in the $-mass region which is due to the meson induced secondary collisions 
channels KK —>■ $ — ► e + e~ and np —>■ $ — ► e + e~. 

In Fig. H] we show the weighted anisotropy coefficients < Bi(M) > for p + Be and Au + 
Au collisions at 10 GeV/ A using the same notations as in Fig.g. The main contribution 
arises from the //-Dalitz decays whereas the contribution from the tu-Dalitz decay is quite 
small due to its reduced cross section. The weighted anisotropy coefficient from pion 
annihilation for p + Be is zero for the same reason. However, the < B^+^~ (M) > for 
Au + Au is zero because the angular distribution of the annihilating pions relative to the 
q AB direction becomes practically isotropic at this energy which leads automatically to 
an isotropic dilepton spectrum. Obviously, this effect is related to our choice of the z-axis 
along the q AB . One can show, for example, that in the Gottfried- Jackson frame, the 
distribution of annihilating pions is anisotropic and < B^+ v - (M) > is small and negative 
at this energy. 

Before presenting our results for SPS energies we like to point out that the transport 
approach HSD provides a good reproduction of the pion rapidity and transverse energy 
distribution at SIS/BEVALAC and AGS energies [13 which ensures that the actual pion 
phase-space densities and thus the 7r + 7r~ annihilation rate should be realistic. However, a 
definite proof for SPS is still lacking, since only rapidity distributions have been presented 
so far in comparison to the experimental data for S + Au at 200 GeV/A in |0|. We thus 
show in Fig. [| the calculated invariant cross section for ir° mesons in the rapidity interval 
2.1 < y < 2.9 as a function of the transverse mass of the pions for S + Au at 200 GeV/A 
in comparison with the experimental data of the WA80 collaboration PB[ . As can be seen 
from Fig. ^ the results of our computation (open circles) are in a good agreement with 
the experimental data (full circles) also with respect to the transverse mass distribution. 

The dilepton invariant mass spectra for p+Be at 450 GeV/A and S+Au at 200 GeV/A 
have been presented in |0|. For 'free' meson masses and formfactors the experimental 
cross section is slightly underestimated for invariant masses 0.3 < M < 0.45 GeV, which 
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may be due to a dropping of the p-meson mass at finite baryon density |T5|, p7[| . In the 
present work we do not explore this idea in more detail and calculate the dilepton spectra 
without any medium modifications. 

In Fig. [| we show the differential dilepton multiplicities for Au + Au at 160 GeV/A, a 
system that has recently been investigated experimentally at the SPS. Again the dominant 
contributions come from the 77 Dalitz decay at small invariant masses and from the pion 
annihilation and direct decay of vector mesons at M ~ m p . 

In Fig. P we present the weighted anisotropy coefficients < Bi(M) > for p + Be 
collisions at 450 GeV/A, S + Au at 200 GeV/A and Au + Au at 160 GeV/A. For p+Be the 
situation is similar to the previous cases at lower energies, the main contribution at small 
invariant masses comes from the rj Dalitz decay; < B n + n -(M) > is approximately zero 
due to the small cross section. For nucleus-nucleus collisions, however, the contribution 
of the pion annihilation channel becomes more essential. As seen from Fig. |6|, contrary to 
the lower bombarding energies, the weighted anisotropy coefficient for pion annihilation 
is negative; an effect related to some differences in the pion dynamics at low and high 
energy as we will discuss below. 



In order to explore this effect we use the result from |T3| that the dilepton anisotropy 
coefficient < B^+^-i^M) > is directly connected with the angular distribution of annihi- 
lating pions W(M, cos 9 n ) relative to q AB : 

b ab (M) _ Jdcos6 n A(6 n ) By%) W(M,cos6„) 
W+7T f d cos 6^ A{6 n ) W{M, cos 6 n ) 

The quantities A{6. K ), B(8 n ) result from a Lorentz transformation from the angle 9 a b to 
Bab and read 

r>ab r)ab 

A(e v ) = 1 + -*±^ sin 2 6 n , B%) = -f^-(3 cos 2 6 W -1). (14) 
1 ^A{6 7T ) 

Here B^ b +w - = — 1, while 6 n is the angle between the pion momentum fi* a in the c.m.s. of 
a + b and the direction of the vector v AB . In the A + B c.m.s. 8 n can be written as 

cos6 w = i , (15) 

\j[M 2 - Ami) (E 2 ~ M2 ) 

where AE = E a — E^, E = E a + are the relative and total energies of the annihilating 
pions with energy E a , E^ in the c.m.s. of A + B. 

In Eq. (|I1J) the 7r + 7r~ angular distribution W(M, cos8 n ) describes the energy and 
angular distribution of the pions that annihilate. Because of simple kinematical relations 
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the distribution of the annihilating pions in AE and E becomes broader with increasing 
bombarding energy. In line with Eq. (|l5|) this leads to a transition from a distribution 
W(M, cos9 n ) peaked at 90° at low energies to a much flatter one at high energies or to a 
positive B coefficient at small energies and a negative one at high energies. 

To quantify our above arguments for the coefficient < B^^- (M) > we show in Fig. [7| 
the distribution W(M,cos9 n ), a contour plot of cos 9 n in the plane (AE,E) at M = 
0.7 GeV (see Eq. ([15])), as well as contour plots of the density distribution W(AE, E) 
for Ca + Ca at 2 GeV/ A, Au + Au at 10 GeV/ A and S + Au at 200 GeV/ A. The value 
of cos 6^ is zero at AE = and goes to ±1 when increasing the absolute value of AE 
according to Eq. First of all one can see a pronounced change in the shape of 

W(M, cos 9 W ) from Ca + Ca at 2 GeV/A to isotropy for Au + Au at 10 GeV/A and to a 
forward - backward peaked distribution for S + Au at 200 GeV/A. In the (£", Ai?)-plane 
the maximum of W(AE, E) is located near AE = at low energy which corresponds to 
cos 9 n = 0. Thus, at low energy the probability for pions to annihilate at angle 9 n = 90° 
relative to the q AB is larger than at 0°, 180°; correspondingly W(M, cos 9 n ) has a maximum 
at 90° (r.h.s). Furthermore, as seen from Fig. |j] the absolute value of the W max (AE , E) at 
AE = decreases with increasing energy and the probability for annihilating pions with 
9 n = 0°, 180° increases. In other words, at high energies the annihilating pions have larger 
longitudinal than transverse momentum in the c.m.s. of leptons relative to q AB (r.h.s.), 
and according to Eq. (|H|) the anisotropy coefficient B A + % - becomes negative. 



4 Summary 

In this work we have presented a fully dynamical study of dilepton production in proton- 
nucleus and nucleus- nucleus collisions from SIS/BEVALAC to SPS energies. Due to the 
parallel ensemble algorithm employed we have, furthermore, calculated dilepton angular 
anisotropies on an event by event basis for the channels rj — > / ye + e~, uj — ► 7r°e + e~, 
A — > Ne + e~ as well as dilepton bremsstrahlung and 7r + 7r~ annihilation. In all reactions 
we find the anisotropy coefficient < B > to be dominated by the 77-decay for invariant 
masses M < 0.4 GeV, while < B > is practically zero for M > 0.5 in p+Be reactions at all 
bombarding energies. This situation changes significantly for nucleus-nucleus collisions 
showing a positive < B > from 7r + 7r~ annihilation of about 0.1 at 2 GeV/A, turning 
approximately zero at 10 GeV/A and becoming negative < B >m —0.1 at SPS energies. 
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We could demonstrate that this energy variation is due to the different pion annihilation 
dynamics. 

Thus, the calculated anisotropy coefficients for nucleus-nucleus collisions support our 
suggestion in Refs. Jl(| [TT], [12| that the dilepton decay anisotropy may serve as an addi- 
tional observable to decompose the dilepton spectra into the various sources and to obtain 
additional information on the reaction dynamics. We note, however, that our transport 
calculations correspond to the full acceptance for the lepton pairs. Especially cuts in the 
transverse momenta affect the anisotropies sensitively. This should be born in mind when 
analyzing experimental data with a limited acceptance. 
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Figure captions 



Figure 1: The weighted anisotropy coefficients < Bi(M) > for p + Be collisions at 
the bombarding energy of 2.1 GeV and Ca + Ca collisions at the bombarding energy 
of 2.0 GeV/A. The "77" denotes the contribution of the //-channel, the "A" labels the 
contribution of the A Dalitz decay, ll pn n the proton-neutron bremsstrahlung, and "7r + 7r~" 
the pion annihilation channel. The solid curves (denoted by "a//" ) show the sum of all 
sources. 

Figure 2: The differential multiplicity cin e + e - / dM for p+Be and central Au+Au collisions 
at a bombarding energy of 10 GeV/A. The ll u — ► 7re + e~" is the uj — ► 7re + e~ Dalitz decay, 
"p", "a/', are the direct decays of the vector mesons. The other notations are the 
same as in Fig. [I]. 

Figure 3: The weighted anisotropy coefficients < B^{M) > for p+Be and central Au + Au 
collisions at 10 GeV/A. The notation is the same as in Fig. [2| 

Figure 4: The calculated invariant cross section (open circles) for tt° mesons in the ra- 
pidity interval 2.1 < y < 2.9 as a function of the transverse pion mass for S + Au at 200 
GeV/A in comparison with the experimental data (full circles) of the WA80 collabora- 
tion m. 

Figure 5: The differential multiplicity dn e + e - / dM for a central Au + Au collision at 160 
GeV/A. The notation is the same as in Fig. |2|. 

Figure 6: The weighted anisotropy coefficients < Bi(M) > for p + Be collisions at 450 
GeV/A, S + Au at 200 GeV/A and Au + Au at 160 GeV/ A. The notation is the same as 
in Fig. 0. 
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Figure 7: (l.h.s.): The contour plot of cos9 n in the plane (AE, E) at M = 0.7 GeV (see 
Eq. (|15])), as well as the countor plots of the density distribution W(AE, E) for Ca + Ca 
at 2 GeV/A, Au + Au at 10 GeV/ A and 5 + Am at 200 GeV/A. (r.h.s.) The angular 
distribution W(M, cos6) for 0.65 < M < 0.75 GeV for the same systems. The solid lines 
represent fits to the calculated distributions. 
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